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FUNCTION OF A SIMPLE WING CONFIGURATION 
IN INCOMPRESSIBLE FLOW 
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SUMMARY 


Measurements of sinusoidally oscillating downwash and lift produced 
on a simple rigid airfoil have been made. Comparison with theory indi- 
cates agreement in order of magnitude and trend, although the scatter of 
results precludes their use for direct verification of Sears' gust func- 
tion. Necessary improvements in experimental techniques are discussed. 

Measurements of statistically stationary random downwash and corre- 
sponding lift on a simple rigid airfoil have been made and values of the 
transfer function between their power spectra determined. An approximate 
analysis has been carried out and the results are compared with the present 
experimental results as well as with those of another investigation. Rea- 
sonable agreement is shown at high reduced frequencies . A more extensive 
measurement program is recommended. 

INTRODUCTION 


Until recently, the emphasis in investigations of aircraft gust loads 
has been on the response to a single gust. Flight test programs, beginning 
with Doolittle's work in 1924 (ref. l) , have been carried out in order to 
determine a reasonable "maximum" gust velocity for airworthiness require- 
ments (refs. 2 and 3). Theoretical considerations, following the funda- 
mental papers of Rhode and Lundquist (ref. 4 ) and of KUssner (ref. 5 ), 
have dealt primarily with methods for computing the transient response 
of an aircraft to such a single gust. 

In reality, however, atmospheric turbulence does not generally con- 
sist of single gusts but is a continuous random process. An airplane 
flying through disturbed air experiences a load history which is a ran- 
domly oscillating function of time. It is therefore apparent that an 
analysis of gust loads based on the concept of a random process instead 
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of that of a single gust would be much closer to the physical nature of 
the problem and could also be expected to yield more realistic and sig- 
nificant results. Unfortunately, the necessary mathematical techniques 
had been lacking until Wiener's generalized harmonic analysis (e.g., 
ref. 6) was applied to practical problems in engineering and research 
(ref. 7). 

By treating the atmospheric turbulence as a statistically stationary 
process and by making use of the response of an aircraft to a purely sinus- 
oidal disturbance, the generalized harmonic analysis provides relations 
between the power spectral densities of the fluctuations of the vertical 
velocity component and of the lift experienced by the airplane. In gen- 
eral, if these power spectra are denoted by $> v and respectively, 

and the sinusoidal transfer function of the aircraft, by r(tu) , and if 
the effects of spanwise variation of the instantaneous downwash are 
neglected, one obtains 

*LL>) = I r ^| 2 ®w^ 

A similar result was utilized by Clementson in 1950 -(ref. 8) for 
deducing power spectra of atmospheric turbulence from the measurements 
of pitch response of aircraft in flight. At the same time, Liepmann 
computed the response of an airfoil to turbulent velocity fluctuations 
(refs. 9 and 10). Since these reports, several papers have been published^ 
about the gust-load problem (refs. 11 to 17) ns well as about the inverse 
problem of atmospheric turbulence measurements (refs . 18 and 19) • 

The present research program at the Massachusetts Institute of 
Technology was motivated by the lack of experimental proof of the valid- 
ity of the new method. Except for concurrent work (ref. 13), there were 
no investigations which could provide such a proof . In fact, even the 
two-dimensional sinusoidal response function, which was computed by Sears 
in 19^1 (ref. 20), had never been satisfactorily verified experimentally. 
It was therefore decided to undertake the following two-phase program: 

(1) To create a sinusoidal cross -velocity fluctuation in a low-speed 
wind t unn el and to measure its characteristics, as well as the character- 
istics of the lift fluctuation on an airfoil of simple geometric configu- 
ration exposed to the velocity fluctuations . From the measurements one 
can then obtain experimental values of the transfer function which are 
directly comparable with Sears' gust function. 

(2) To create a random but statistically stationary velocity dis- 
turbance field in a low-speed wind tunnel and to measure the power spectra 
of fluctuations of velocity and of lift on a rigid airfoil. From the 
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measurements one can again obtain experimental values of the amplitude 
gain of the transfer function, and a comparison with theoretical results 
can be made. 

During the course of the investigation, several problems which were 
secondary to the main objective turned out to be of importance. One of 
these is the effect of spanwise variation of the turbulent velocity flue 
tuation, which modifies the sinusoidal transfer function much more than 
was originally expected. Another problem to which considerable effort 
had to be devoted is that of the measurement of statistical quantities 
with sufficient accuracy. 

This research was conducted under the sponsorship and with the 
financial assistance of the National Advisory Committee for Aeronautics. 


SYMBOLS 


area of force-sensitive wing section, sq in. 


aw* 

Ry — Rg 

hot-wire overheating ratio, 

R 

a(a>,o^) 

filter response curve 

b 

span of force-sensitive wing section 

Ci 

0 00 

/v / cos 

cosine integral, Ci(,|) - - J ^ 

c L 

lift coefficient 

c 

chord, in. 

F(flCb) 

filter factor, radians/sec 

f 

filter output, mv 

g 

span-effect correction factor 

OJ 

•s 

c\T 

£ 

Hahkel functions 

h 

_2 

lift influence function, in. 
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h b 
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Io 

J o> J l 

Ko,K! 

k 

k l 

k 2 

L 

l 

Q 

q 

R e 

R* 

r 

S 

Si 

s 

T 

T c 


wing oscillation amplitude, in. 

strip-theory lift influence function, in .“ 2 

_p 

total lift influence function, in. 

hot-wire current, ma 

hot-wire calibration constant, ma 

Bessel functions of first kind 
modified Bessel functions of second kind 
reduced frequency 

streamwise wave number, radians /in. 
spanwise wave number, radians /in. 

lift, lb 

streamwise distance, in. 
filter-band-width parameter in 

dynamic pressure, psi 
hot-wire cold resistance, ohms 
hot-wire working resistance, ohms 

hot-wire output tapping ratio 
Sears 1 gust function 

sine integral, Si(|) = f - &*> 1 

Jq ^ 

thermocouple output, rnv 
integration time, sec 
thermocouple time constant, sec 


|a(o),uo) 


1 + 

V , 
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t 

time, sec 

u 

stream velocity, mph and in. /sec 

u 

streamwise velocity fluctuation, in. /sec 

V 

crosswise velocity fluctuation, in. /sec 

X 

streamwise coordinate, in. 

y 

spanwise coordinate, in. 

a 

angle of attack, radians 

|3,e 

dimensionless variable in g 

r 

transfer function in one variable 

? 

transfer function in two wave numbers 

7 

dimensionless thermocouple time constant 

T] 

spanwise distance variable 

A 

turbulence scale, in. 

V 

hot-wire output voltage, rav 

1 

streamwise distance variable, in. 

P 

dimensionless integration time 

a 

root-mean-square error 

T 

time variable, sec 

$LL 

lift power spectrum 

$ 

s l s l 

power spectral density of input 

^w 

power spectral density of v in one variable, 
( in . /sec ) 2 (radians /in . ) ~ 1 

$ 

w 

power spectral density of v in two wave numbers, 
(in. /sec) 2 (radians /in. )" 2 

<Pll 

lift correlation function 

T’ss 

normalized correlation function of s(t) 
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cp^ correlation function of v in one variable, (in. /sec) ^ 

correlation function of v in two variables, (in. /sec) 2 

\(r hot-wire angle, radians 

ft reduced circulation amplitude, in. /sec 

cd circular frequency, radians/sec 


GENERAL EXPERIMENTAL EQUIPMENT 


Wind Tunnel 

All measurements were performed in the 5~ by 7^ ~ foot low-speed wind 

tunnel of the M.I.T. Aeroelastic and Structures Research Laboratory. This 
tunnel has a maximum speed of approximately 100 mph. However, in order to 
avoid excessive overheating of the tunnel during long runs most of the 
present measurements were made at a speed of 50 mph. 

Since this wind tunnel was originally built for flutter investiga- 
tions, no special care was taken in its design to obtain a low turbulence 
level. As shown in figure 1, both the contraction ratio and the distance 
from the turning vanes to the test section are relatively small. A tur- 
bulence level of the order of one-quarter of 1 percent root mean square 
of the free-stream velocity was measured at a speed of 50 m Ph with an 
uncompensated 0.00015-inch tungsten hot-wire, which has a frequency- 
response cutoff around 200 cps. 

The tunnel velocity was generally measured by means of a pitot-static 
probe and an inclined-tube alcohol manometer. In hot-wire calibrations, 
very low velocities (of the order of 3 fps) were needed and the vortex- 
frequency technique described in reference 21 was successfully used. In 
the overlapping region of the ranges of the two velocity measuring devices 
an agreement within about 2 percent was generally obtained. 

Hot-wire surveys of the flow indicated satisfactorily uniform condi- 
tions in the region where measurement apparatus was located. 


Measurement of Velocity Fluctuations 

As will be described in the succeeding sections, additional velocity 
fluctuations were introduced into the mean flow by mechanical devices. 

The component of these fluctuations normal to the test airfoil was meas- 
ured by means of a directionally sensitive hot-wire probe. This probe 
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consisted of two approximately 0.05-inch-long and 0.00015-inch-diameter 
tungsten wires in the form of an x. The probe was mounted on a vertical 
stand and could be both moved vertically and turned in the horizontal 
plane by remote control. The horizontal position of the stand could be 
adjusted between runs. A special rigid-base structure attached directly 
to the building walls was brought in through holes in the tunnel floor 
and ceiling in order to minimize vibrations of the probe. Figure 2 shows 
the probe in position. 

Facilities for manufacture and repair of the hot-wire probes were 
set up at the laboratory. After some initial difficulties, no trouble 
beyond occasional breakage and dirt collection was experienced. 


Electrical Circuit 

It is difficult to produce two exactly identical wires for a double- 
wire probe. Consequently, in order to obtain pure lateral-velocity 
response, the heating currents in the two wires must be individually 
adjustable. Therefore, a circuit consisting of two separate Wheatstone 
bridges was built (fig. 3). Currents were measured by reading the volt- 
age drop across a precision resistor on a Rubicon potentiometer. Bridge 
balance was observed on a Weston sensitive galvanometer. 


Calibration of Hot-Wire Probes 

The conditions for operating a double oblique -wire probe with pure 
cross-velocity sensitivity are derived in appendix A. First, an individ- 
ual calibration of each wire is needed for determination of the calibra- 
tion constants, from which the current and voltage divider settings can 
be calculated. Such calibrations were made normally at four velocities 
and two values of the overheating ratio. Calibration of the cross -velocity 
sensitivity was obtained by turning the probe in a horizontal plane through 
known angles . Independence of fluctuations in the streamwise direction was 
finally ascertained by observing the probe response while tunnel speed was 
varied within ±10 percent about the operating condition. Typical single- 
wire and cross -velocity calibrations are presented in figures 1+ and 5> 
respectively. 

Satisfactory linearity and repeatability were generally obtained in 
the calibrations. Occasionally it was found, however, that after a long 
time of operation the calibration changed slightly because of collection 
of dirt on the wires. After cleaning the probe, the original calibration 
was usually obtained, in some cases with remarkable accuracy. The wires 
were cleaned with ether every day before runs were started, and it is not 
believed that serious errors could have been introduced into the measure- 
ments by the dirt deposits. 
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Since the frequency range of interest in the present investigation 
extended only up to about 150 cps, no compensating circuits were used 
for the 0.00015-inch tungsten wires. 


Measurement of Lift Fluctuations 

Since it was desired to obtain measurements of lift produced by tur- 
bulence on a simple airfoil configuration, a constant-cross-section wing 
spanning the 5-foot tunnel height was installed on the same rigid 
supporting structure as the hot-wire stand. The wing chord was chosen 
to be 1 foot in order to provide a reasonably high Reynolds number 
(500,000 at 50 mph tunnel speed) . The profile is NACA 0010. 

The main structure of the wing consisted of two full-span 1- by 

2 . — inch steel beams and two i - by 2^ -inch steel flanges, bolted 
16 8 4 

together. Profile shape was obtained by attaching balsa blocks to the 
steel parts. In the center of the wing, a 4-inch-span section was cut 
off, except for the two main steel beams, and replaced by a load-sensitive 
element mounted on strain-gage pickups. Figure 6 shows details of the 
wing. 


The load-sensitive section consisted of two hollowed balsa blocks 
joined together by thin curved magnesium plates. One of the plates was 
supported by three strain-gage load cells in such a manner that the sec- 
tion "floated" about the throughgoing steel structure. Thus, a lift load 
on the section was indicated as the sum of the signals from the three 
pickups. The load cells were manufactured by the Dynamic Instrument Co. 
Each one consisted of unbonded strain gages, arranged as a four-active- 
arm Wheatstone bridge, supporting a steel pin which protruded through 
the end of the cylindrical-shape pickup. Dimensions of the pickup were 
small enough to permit complete internal mounting within the 1.2-inch- 
thick airfoil (see fig. 6). Connection to the magnesium plate was made 
by means of a brass bushing and soft soldering. Current to each load 
pickup was supplied by a separate dry cell, and the outputs were con- 
nected in series. Figure 7 shows the circuit diagram and figure 8, a 
typical calibration of one pickup. Before each series of runs, a cali- 
bration of the complete system was made and its independence of the point 
of application of the load was checked. 


The principal consideration in the design of the lift sensing sys- 
tem was to obtain a mechanism rigid enough to avoid excessive displace- 
ments and resonance peaks in the important frequency range (below 200 cps ) . 
Both of these purposes were satisfactorily accomplished; the displacements 
were of the order of 0.001 inch or less, and a natural frequency of 330 cps 
was obtained. The dynamic response was tested by means of an electromag- 
netic shaker. The floating wing section and the pickups were mounted on 
an exact duplicate of the local wing-span structure (shown in fig. 6) and 
shaken at constant acceleration throughout the desired frequency range. 
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In the first tests, a strong resonance peak was obtained at around 200 cps. 
It disappeared completely, however, after stiffening the structure of the 
floating section by cementing strips of magnesium sheet to the sides of 
the rear balsa block. 

A serious problem was presented by the vibration of the entire wing 
structure due to the turbulence acting on it. The wing had a natural 
bending frequency of only approximately 25 cps and considerable vibration 
was therefore present at the load cell mountings . Instead of trying to 
stiffen the wing structure by guy wires or otherwise, it was decided to 
install a separate pickup inside the floating section so that it would 
respond to only the structural vibrations. By suitable adjustments of 
the feeding voltage and of a dummy mass attached to the pickup, the 
dynamic responses of the two systems could be made nearly identical. 

Thus, if the signal of the dummy system was reversed and connected in 
series with the main signal, the structural vibrations would cancel out 
and a signal indicating the pure aerodynamic force could be obtained. 

The electromagnetic shaker was again used for verification of the cancel- 
lation performance. The measured dynamic response curves, referred to 
the uncompensated value at zero frequency, are presented in figure 9* 

It should be noted that, while the cancellation is very effective at 
low frequencies, considerable errors appear at certain frequencies and 
especially around the resonance. Since, however, there is an appreciable 
reduction in the support vibration response throughout the frequency range 
under investigation (even in the worst case at 150 cps the cancellation is 
about 75 percent), it was not considered necessary to improve the perform- 
ance by closer matching of damping characteristics or other methods. 

As far as support vibrations in pitch and roll are concerned, their 
effects can theoretically be shown to cancel out if the center of gravity 
of the floating section is placed in the exact geometric center of the 
triangle formed by the three pickups and if the dummy pickup is located at 
that same point. The first condition was satisfied exactly by use of small 
lead weights in the front balsa block. In the second one, a small error 
was necessary because of space requirements. In any case, in view of the 
relatively high torsional rigidity of the wing, the pitch and roll vibra- 
tions are believed to be of secondary importance. 

The static aerodynamic response of the lift sensing system was 
obtained experimentally by turning the wing about a vertical axis to 
known values of angle of attack. The resulting lift-curve slope was 
0.08 per degree with excellent linearity through the ±7° range covered 
in the measurement. 


Auxiliary Equipment 

A two -channel amplifier was designed and built for use with both 
hot-wire and lift signals. The amplifier operates on the push-pull prin- 
ciple and has a one-side maximum gain of 2,b00, with a noise level of 
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approximately 100 microvolt referred to the input. The noise was, how- 
ever, found to consist largely of 120-cps pickup and could he avoided in 
the spectrum measurements . Frequency response of the amplifier was flat 
from 0.5 cps to several kilocycles. Circuit diagram, amplitude calibra- 
tion, and frequency response are shown in figures 10, 11, and 12, 
respectively. 

By means of a switching circuit, amplifier inputs could be selected 
from either external sources or from the bridge circuit, such as the hot- 
wire signals, their difference and sum, and so forth. The amplifier out- 
put was available as push-pull or low-impedance cathode follower signals 
or could be fed into a heater -thermocouple unit through a driver stage 

(fig. 13). 

The heater -thermocouple unit produces a smoothed square of the input 
signal. Calibrations for steady direct current or sinusoidal signals are 
given in figure 14- . For random signals the output fluctuates, making 
visual reading of the mean-square value difficult. To minimize this 
fluctuation, units with a time constant of 5 seconds were obtained from 
the Cambridge Instrument Co. in England. Nevertheless, it turned out to 
be necessary to average the thermocouple output more accurately by inte- 
grating it over a long period of time (see section on turbulence measure- 
ments) by means of a Miller type integrator. The integrator consisted of 
a Reeves direct -current amplifier and a 1-microfarad condenser connected 
between input and output; it has the calibration shown in figure 15 • 


RESPONSE OF RIGID AIRFOIL TO SINUSOIDAL GUSTS 


Experimental Techniques 

The sinusoidal velocity fluctuations were generated by a 1-foot-chord, 
5-foot-span wing which was rigidly mounted on a horizontally oscillating 
steel rod at the upstream end of the test section (see fig. 16). The maxi- 
mum amplitude of the translatory oscillation was approximately 2 inches and 
the practical frequency range extended up to 10 cps, although at the higher 
frequencies the amplitude had to be considerably decreased owing to con- 
siderations of wing strength and rigidity of the oscillator mechanism. The 
wing had an NACA 0012 airfoil and a magnesium-skin balsa structure . 

The velocity fluctuations and lift produced on the stationary wing 
were measured by the hot-wire and strain-gage equipment described in the 
preceding section. The frequency of the moving-wing oscillation was 
measured by means of a microswitch and a Hewlett-Packard electronic 
counter . 

The root-mean-square sinusoidal gust velocity was found to be less 
than one-half of 1 percent of the free-stream velocity. Although both 
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the velocity and lift signals were well within the useful range of their 
respective sensing units, they could not he satisfactorily distinguished 
from the effects of the relatively high free -stream turbulence level of 
the wind tunnel. It was therefore decided to filter the signals elec- 
trically by means of a General Radio 7^2-B vibration analyzer . This 
instrument has a narrow-band-pass characteristic, with a value of the 
Q-parameter of about 80. 

Use of the narrow-band filter greatly reduced the effects of back- 
ground turbulence. However, a new difficulty appeared because of a drift 
in the oscillator frequency of the order of ±0.5 percent of its mean value. 
With the narrow-band -pass filter, this drift caused fluctuations in the 
measured signal of about ±30 percent. It was estimated that in order to 
decrease this figure to a satisfactory level a frequency drift of only 
±0.1 percent would be permissible. The obvious remedies would be either 
to obtain a better frequency control for the oscillator or to use a fil- 
ter with a wider flattop characteristic curve. Both of these measures 
would, however, have required design and construction of precision machin- 
ery or accurate electronic circuits and involved considerable effort and 
expense beyond the resources of the present program. 

It was originally planned to obtain measurements of the phase angle 
as well as of the amplitude gain of the Sears' function. Unfortunately, 
the phase shift in the band-pass filter unit turned out to depend strongly 
on the signal amplitude and accuracy of timing, thus making the phase 
measurement impractical. After the experiences described above, it was 
concluded that the only way to obtain satisfactory measurements of the 
Sear's function would be to perform the experiments in a tunnel with a 
low background— turbulence level and with an oscillator mechanism which 
would be capable of producing large downwash amplitudes. It is believed 
that significant improvements can be made in both of these items in the 
present wind tunnel, but only after considerable effort and expense. 

A question was raised regarding the possibility that even with 
absolutely steady oscillator frequency the sinusoidal downwash pattern 
downstream may have a tendency to "contract" and "expand" in response, 
to low-frequency disturbances in the streamwise velocity, thus producing 
hot-wire and lift signals which would exhibit a frequency drift. The. 
present instrumentation could not provide conclusive results about this 
suspected phenomenon; it is, however, believed to be worthy of further 
investigation. 

In spite of these difficulties, it was decided to make such meas- 
urements as were possible in order to gain experience for further work. 

The measurements for downwash and lift amplitudes were made in separate 
runs with the hot-wire mount and the rigid wing in approximately the 
same position. Actually, because of structural reasons, the hot-wire 
was located about l/2 chord upstream of the wing leading edge; however, 
since phase measurements were not attempted, the resulting error is 
negligible . 
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Measurement of Sinusoidal Velocity Fluctuations 

After filtering, the hot-wire signal was amplified and its mean 
square was determined by means of the heater-thermocouple unit. Since 
considerable fluctuation was observed in the thermocouple output, inte- 
gration over several minutes was necessary. The system diagram is pre- 
sented in figure 17. 

By shutting off the oscillator system, it was possible to deter- 
mine the mean-square value of the combination of remaining background 
turbulence and amplifier noise. Subtraction of this value from the mean 
square of the total signal should then provide the mean square of the 
pure sinusoid. This correction cannot, however, account for any effects 
of the frequency drift. Figure 18 shows sample recordings of the hot- 
wire signal taken by means of a Brush oscillograph. 

For calibration purposes, a known sinusoidal signal was injected 
into the filter input. Frequency of the calibration signal was deter- 
mined by comparing it against the oscillating -wing frequency on a 
Hewlett-Packard counter. Considerable scatter is shown in the results 

(fig. 19). 


Measurement of Sinusoidal Lift Fluctuations 


A somewhat different technique was used for analysis of the lift 
signals. For reasons which are not completely understood, a consider- 
ably lower noise level was observed in the lift signal than in the hot- 
wire output. In part this is believed to be caused by some high- 
frequency noise leaking through the filter, as was observed to happen 
to some extent in the hot-wire measurements. Because of the aerodynamic 
response characteristics of the wing, much less such noise is present in 
the strain-gage signal. 


The amplitude of the lift signal could now be obtained as an average 
of several visual readings on an oscilloscope at short intervals (fig- ) 
The calibration signal, its frequency determined as in the case of the 
hot-wire measurements, was also observed on the oscilloscope. 

The results (fig. 21) show somewhat less scatter than the downwash 
measurements . 


Theoretical Considerations 

According to Ashley (ref. 22), the downwash in the plane of the wake 
behind a two-dimensional wing in sinusoidal oscillation can be expressed 
as follows, provided that the point of observation is sufficiently far 
from the wing (distance l from midchord) : 
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v(l,t) = - 


2rt 


( e -lk 

■L j(2Z/c) - 1 


\(2l/c) + 1 



where Si and Ci are the sine and cosine integrals , respectively , and 
the reduced circulation amplitude is 


ft = 


Li 


1 + f ' 

-f7V 0 (S')«* 


jtik 


H^ 2) (k) + iH^ 2) (k) 


and where v 0 (|') is the local amplitude of the wing oscillation velocity 
as a function of chordwise position | ' = x In the present case of 

translatory motion with amplitude h. 


Vq(I') = -icoh 


and 


a = 


to ( s /t) 

4 2) (k) + 1 Hq 2) (k) 


By substitution, the complex downwash amplitude becomes 


v(Z) = 



Si - 1 



Ik 
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A plot of the absolute amplitude of the downwash angle 

in figure 19 for the present case of 2l/c = 13* Since 
curve slope of the wing is only 4-58 per radian instead 
2rt per radian, a corresponding correction has been made 
curve. 


(5 )/(tt) is s iven 

the measured lift- 
of the theoretical 
to the plotted 


The sinusoidal lift produced by this downwash on another two- 
dimensional wing of the same chord has an amplitude (assuming that the 
downwash over the wing is constant and equal to the calculated value at 

midchord) of 


L = 


v(D 

U 


(2*) [s(k)] (qA) 


where S(k) is the Sears' gust function. 


The absolute dimensionless amplitude 


L 1 _ 

cl 

$ qA| 

2h/c 


is plotted 


in figure 21. Since the measured steady-state lift-curve slope of the 
lift-sensing wing is dC^/da =4-58 per radian, this value was again 
substituted for the theoretical 2jt per radian. 


Discussion 

Because of the scatter of the measurements, it was not possible to 
compare the results directly with the results of Sears' gust function. 

The measured lift, however, agrees reasonably well with the theoretical 

prediction. 

The equipment developed for measuring downwash and lift was found 
to be satisfactory. However, a much lower ratio of background turbulence 
to sinusoidal downwash in the wind tunnel and a better frequency control 
of the oscillator are required for successful completion of conclusive 
measurements . 


RESPONSE OE A RIGID AIRFOIL TO RANDOM GUSTS 


Experimental Technique 


a 


For measuring the random 
rigid airfoil, the equipment 


downwash and the associated lift response of. 
described in the preceding sections was again 
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used. Both the hot-wire and the lift-sensing wing were placed in the 
center of the test section approximately 2 feet apart in order to permit 
simultaneous measurements of lift and downwash. This arrangement elimi- 
nates errors due to differences in operating conditions in separate runs. 
It is not believed that interference effects large enough to offset this 
gain could have been introduced into the turbulence field. 

The downwash fluctuations were generated by a coarse grid installed 
approximately 10 feet upstream of the measuring units (fig. 22). The 
mesh size of the grid was 6 inches, which was considered to be a maximum 
for maintaining a reasonably homogeneous turbulence field in the test 
area. The grid members were 1 inch wide with a square cross section. 

Unfortunately , the turbulence scale obtainable with this grid is, 
at most, of the order of the mesh size, that is, less than a half -chord 
of the wing. In atmospheric turbulence, scales much larger than the 
chord length of an average airplane are generally encountered (ref. 8). 
Naturally, a closer matching could be obtained by decreasing the model 
size; this, however, would result in very low Reynolds numbers, with 
associated adverse aerodynamic effects, as well as in instrumentation 
difficulties . 

A hot-wire survey of the turbulence level in the test area showed 
homogeneity of the root -mean-square fluctuation to tunnel velocity ratio 
within approximately 5 percent. The average value of this ratio was 
approximately 0.03* A background-turbulence level of about one-fourth 
of 1 percent root mean square was observed in the empty tunnel. The 
turbulent fluctuations are therefore actually generated by two separate 
sources, although the contribution of the grid is dominating. Since the 
power spectral density of the total turbulent downwash is measured, the 
dual origin does not present serious complications. 


Recording and Analysis of Data 

In order to avoid excessively long running times with associated 
changes in tunnel operating conditions, both the hot-wire and lift sig- 
nals were recorded simultaneously on magnetic tape by means of an 
Ampex 306 two-channel recorder. The recording time was generally sev- 
eral minutes and was determined on the basis of accuracy considerations 
which will be discussed subsequently. The system diagram is presented 
in figure 23* 

The power spectra were determined by playing the recorded signal 
into the thermocouple squaring unit through a narrow-band -pass filter set 
at different frequencies. A General Radio 75^-B vibration analyzer was 
used as the filter. This instrument has a Wien bridge circuit with a 
band -pass parameter Q = 80 and a range from 2.5 to 750 cps. 
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Some difficulty was experienced in the calibration of the vibration 
analyzer. If the frequency -response curve of the filter at central fre- 
quency o^j is a(a>,nt,), the output signal sp is related to the power 
spectrum of the input sp as follows : 


r 00 2 

= J |a(a>,ab)| ® Sl s 1 ( ai ) dal 


If $ is constant or, in case of a symmetrical (with respect 

s l s l 

to 0^) |a(o),<Xb)| curve, linear over the effective filter range, one 

can simplify 




s l s l 


(cd 0 )F(u» 0 ) 


Thus, the filter response factor F(ot,) is needed for quantitative eval- 
uation of power spectrum by the method described above. A careful meas- 
urement of the response curve |a.(o>,a^- ) ) | was made for selected frequencies 
by actually determining the attenuation of a steady sinusoid in the filter 
range. Then F(<xp) was computed by manual integration of |a(o>,ofc)| . 

Unfortunately, it was found later that the calibration of the vibration 
analyzer was very sensitive to any changes in battery voltages and com- 
ponent characteristics. Since, for instance, batteries and tubes had to 
be replaced during the course of the measurements and a complete recali- 
bration every time would have been impractical, it was decided not to 
attempt accurate absolute measurements of power spectra. 

The main purpose of the experiments, to obtain values of the transfer 
function connecting the input and output power spectra, could be accom- 
plished without knowledge of F(oi 0 ) . Since the hot-wire and strain-gage 

outputs were recorded on magnetic tape, it was possible to run both sig- 
nals through the filter successively, so that changes in the filter 
between the runs were very unlikely. Thus, even though the measured 
mean— square values could not be absolutely evaluated to give the two 
power spectral densities, the unknown proportionality factors were iden- 
tical and canceled out in computation of the transfer function. "Hie 
experimental results presented in figure 2k were determined in this man- 
ner. The physical frequency range of the measurements is from 3 to 
150 cps. An additional dynamic -response correction corresponding to 
figure 9 was applied to the lift measurements at high frequencies. 
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Accuracy of Measurements 

The errors in the present measurements can be divided into the fol- 
lowing two categories : 

(1) Measurement errors, arising from changes in the equipment response 
inaccuracies in calibration, reading errors, and so forth. It is difficult 
to estimate the total probable error of a system as complicated as the one 
used here. The accuracy of voltmeters and other instruments used in cali- 
brating individual parts of the system is generally within ±1 to ±3 percent 
Repeatability of the hot-wire calibrations is of the same order, and inev- 
itable changes occur in the electronic circuits because of tube aging, line 
voltage fluctuations, and so forth. Therefore, a total error of several 
percent has to be considered probable. 

(2) Statistical errors, arising from the fact that the power spectral 
densities are measured as finite time averages. Since the power spectral 
density can be considered an infinite time average of the filtered and 
squared signal, a measurement over a finite time has to suffer in accuracy. 
Fortunately, it is possible to estimate this error on the basis of results 
presented by Rice (ref. 7) • Such an analysis applicable to the present 
measurements is presented in appendix B. 

The calculated root -mean-square error curves in figure 25 were used 
in choosing averaging times at each selected frequency. In general, a 
root-mean-square error of 3 percent was planned for in the power -spectrum 
measurements . The experimental value of the transfer function, computed 
as the ratio of two such spectrum measurements, could then be expected 
to have a statistical root -mean-square error of the order of 5 percent. 

One has to realize, however, that these numbers refer to the root-mean- 
square error of an infinite ensemble of measurements performed under 
identical conditions. All that can be said regarding an individual meas- 
urement is to state its probability of being within specified limits of 
error. For instance, if the probability density distribution of the 
transfer function is assumed to be Gaussian (approximate analysis indi- 
cates that this may not be an unreasonable assumption), a 5-percent root- 
mean-square error implies that a single measurement would have an approx- 
imately 90-percent probability of being within ±10 percent of the correct 
value . 

A series of individual measurements at different frequencies would 
show scatter of the order of the root -mean-square error. In the present 
case, it is not possible to say definitely whether the few points at low 
frequencies reflect an unduly large scatter or an anomalous behavior of 
the transfer function as compared with the theoretical estimate. At 
higher frequencies, where the computed averaging times were also gen- 
erally exceeded for practical reasons, the scatter is of acceptable 
magnitude . 
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The best way to increase the reliability of the present results 
would be to obtain a large number of power -spectrum measurements at each 
selected frequency. The average of such measurements would be essen- 
tially equivalent to a result obtained by using a longer, otherwise 
impractical, averaging time. The statistical accuracy would be corre- 
spondingly improved and some measurement errors of random nature would 
also be minimized. 


The measured characteristics of the turbulent downwash fluctuations 
and the corresponding lift on an airfoil are the power spectral den- 
sities i> w and <6^ in terms of the reduced frequency k = <nc/2U or 

the streamwise wave number kq = ca/U. The transfer function T is then 
obtained from these two spectral densities as follows : 


If the velocity disturbance field were variable in the stream direc- 
tion only, so that the instantaneous angle of attack were constant along 
the span, the transfer function r(kq) would be identical with the trans- 
fer function for a wing passing through a sinusoidal gust field. For a 
wing of constant cross section and infinite span this is the Sears gust 
function (ref. 20): 


Unfortunately, in reality random velocity disturbances in the atmos- 
phere as well as in a wind tunnel possess statistical characteristics 
which are (in most cases) very much the same in any direction the turbu- 
lence field may be traversed. Consequently, in trying to obtain a real- 
istic estimate of the transfer function, one has to account for the span- 
wise variation of the instantaneous angle of attack. This brings in a 
new variable, the spanwise wave number k 2 * An approximate two -wave - 
number analysis will now be carried out for the geometrical configuration 
used in the present experiments. The analysis follows essentially the 
line of thought presented by Liepmann (ref . 15) • It should be emphasized 
that the general procedure is by no means restricted to the simple wing 
geometry and turbulence spectra considered herein. 


Theoretical Considerations 



J 0 (k)Kq(ik) + iJq(k)Ko(ik) 
Kq(ik) + Ko(ik) 
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Using the coordinate systems defined in figure 2 6 and assuming the 
validity of Taylor's hypothesis for the motion of the wing through the 
turbulence field v(x,y), one can write the lift dL on a chordwise 
strip of width dy at y as 

00 

dL(x,y) = -j— — -^j- <ty jj j 7 ! ~ y |) v ( x - 

where [t] - y|) is the influence function for lift per unit span 

at y produced by a unit downwash impulse at position (|,T}) with respect 
to the wing. The total lift as a function of the position of the wing in 
the x-direction is then 


L(*) - ^ 


da 


ob/2 

l-b/2 


00 

dy - y|)v(x - 6,Tj)di d n 


At this point it is convenient to introduce an integrated influence 
function 


rb/2 


hb(l ’” ) ■ti b/2 h ( 5 ’i n - y 0 


<ty 


which represents the lift on the whole lift-sensitive section of span b 
due to a unit impulse downwash at (|,T|) (see fig. 2.6). After substitution 
one obtains 


L(x) = ^ JJ hfrit ,i\)v(x - i,T))d| dtj 


To bring in the power spectral density, one forms the correlation 
of l(x ) , which is taken as an ensemble average over L(x)L(x + Ax): 

L(x)L(x + ^) = bb(6l>ni)bb(8 2 » 7 l 2 )v(x - Sl»^v(x + AX - «2 dT >l dT >2 

—00 
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If "the turbulence field is homogeneous in the xy -plane, the ensemble 
average over the velocity product is independent of x and can be replaced 

by the two-dimensional autocorrelation cp^yjAx + - £2)' (^1 " ^2)] ' Simi- 

larly, the lift correlation can be expressed independently of x, giving 


2 °° 

^(Ax) = fj[[ h b (i 1 ,t 1 )h b (l 2 ,n 2 )$- v |A>c + (ip - l 2 ),(m - d| 2 dn 2 


Now, can be written in terms of its Fourier transform, the 

power spectrum $ v - v .^kp,k2 ^ } as 


<Pw 


. ci rr ~ / \ i [^i( Ax+| i" | 2) +k 2( Ti i _Ti 2) 

Ax + (^l“ ^2) > ^2) = JJ ®w( k : 


dk^ dk 2 


By substitution and interchanging the order of integration 

<Pnfc*> - ^r) 2 /.” elklAX -/I jj' h b(l 1 .1i) exp ifal! + 


cx> 

k 2 T 1 l) d ^l dT ll JJ h b(^ 2 ^ 2 ) ex P _i ( k 1^2 + ^2)^2 dT l 2 


If one now defines a two -wave -number transfer function 


00 

r(ki,k 2 ^ = j^T h- b (|,Ti)e 1 ( kl ^ +k 2 T 1 )d| dr] 


the lift correlation can be written 


dkr 


dkn 


9ll( A x ) = J ^ e lklAx fjtwfc l ,k 2) r ( k l ,k 2) 


'dk 2 


dk n 
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From the uniqueness of the Fourier transform it follows then that the 
lift power spectrum must be 


^LL^l) 



Hence, instead of Sears' function, one obtains a two -wave -number transfer 
function r^k-^kpj and consequently also has to deal with a two -wave - 
number turbulence spectrum $> vv ^k b ,k2^ • 

In trying to obtain an analytical expression the result of which 
could be compared with the experimental result, it is seen from the above 
that only $T,T,( k l) i- s measured directly. The measured turbulence spec- 
trum <& vv (k 1 ) is related to through integration as follows: 


$ w( k i) = $w( k l> k 2) dk 2 


Furthermore , there are no accurate aerodynamic theories for a calcu- 
lation of Consequently, one is forced to make the following 

simplifying assumptions in order to establish a comparison with the 
experimental results : 

(l) It is assumed that the lift influence function follows the aero- 
dynamic strip theory; that is, the lift at any spanwise position y on 
the wing depends only on turbulence at T) = y. This leads to 

h^i, |t) - y|) = h 0 (|) 5 (ri - y) 

h b = ^ (hi = 1) 

h b = 0 (hi >1) 

The strip-theory assumption is certainly not correct, because it elimi- 
nates all effects due to the semi-infinite wing extensions. The resulting 
approximate formulas are therefore identical with those derived with the 
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same assumptions for a finite-span rectangular wing. Too high lift values 
can be expected, especially at a low aspect ratio. However, it is proba- 
ble that in reality the "tails" of the influence function (fig. 26) die 
out rather rapidly, so that the error will decrease with increasing span. 
In any case, it is believed that no fundamental properties of the result 
would be lost because of application of strip theory. 

(2) It is further assumed that the lift response of each strip is 
characterized by the Sears' function. Thus, bo(|) will be the Fourier 
transform of Sears' function, and one can write 




S(k) 


By substitution, the two -wave -number transfer function then becomes 





f h 0 (t)e lkll d| 

/—oo 


r 




The lift power spectrum can now be written 





dk 2 


(3) In order to bring in the measured one -wave -number turbulence 
spectrum ® vv (k 1 ), the following span-correction factor is defined: 




dk.2 
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The lift power spectrum becomes then 



In this expression everything else is known or measured except 
g(b,k^J. Actually, nothing has been gained so far, because calculation 
of g^b,k-^ requires information regarding the two-dimensional turbu- 
lence spectrum ®(k 1 ,k 2 ). However, g(b,k) contains only integrals of 
the spectrum and consequently may not depend too much on the exact expres- 
sion for ® V v^ki,k 2 ) . The possibility of using a simplified expression 
for computation of g^b,k^ becomes apparent; even then, however, some 

connection must be retained to the actual measured spectrum through one 
or more free parameters . 

Experimental results (for instance, refs. 23 and 2 ^) indicate that 
grid turbulence can be reasonably described by a simple spectrum which, 
assuming homogeneity and isotropy, leads to the following form of the 
two -wave -number power spectrum for the fluctuation component normal to 
the main stream (see also ref. 15)2 



The integral scale A, formally the integral over the corresponding nor- 
malized correlation function for streamwise fluctuations, represents a 
characteristic length of the turbulence and is the free parameter to be 
determined by comparison with the measured spectrum. In this matching, 
use can be made of the integral 


r°° / \ / x A V 2 (l + 3 A 2 kl 2 ) 

NoN - LN o( ki ' k2 ) dk 2 - 5 (l + A2ki if 

which can be compared directly with the measured ( k i) • 
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Substitution of 
finally 



into the expression for 



yields 


g(b,ki,A) = ±—r 

v ' 1 + 3P 2 


(1 - 5P 2 ) 


1 - 2€Kp(2e) 


- 2K Q (2e) + ^ f K 0 (2e')de' 

6 Jo 


where Kq and Kp are modified Bessel functions of the second kind and 


e = — \|l + A^k 

2A v 


V 


(3 = Akp 

A graphical presentation of g(b,kp,A) is best given in terms of the 

dimensionless parameters A/b and ^k = Such a plot is presented 

in figure 27 , while figure 2k gives the complete transfer function 
| S(k) | 2 g^b,^,A^ for the present experimental configuration at A/b = 1 

and l/2. It should be noted that even if A/b-^oo, Sears' function is 
obtained as a limiting case for all values of k only when b/c— >0, 
that is, for an infinitely narrow strip. 

It can also be shown that, for small values of b, g(b,kp,A^ can 
be expressed in power series as 


g = 1 - Constant x b 2 + 


(b « A) 


Thus, the power spectral density and its integral, the mean-square lift, 
are proportional to the square of the wing area for small values of span. 
This result is in accordance with the general behavior of linear systems 
responding to turbulent excitation. For high values of b/A, however, 
one obtains the following asymptotic expression: 



1 

1 - 3A 2 kp 2 

1 + 3A 2 kp2 

jSfk + A2kl2 ) 


3rtA 2 kp 2 
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This means that for large-span wings at low wave numbers, the span correc 
tion behaves proportionally to l/b 2 , whereas at high wave numbers it 
becomes approximately l/b. If one then considers the lift power spec- 
trum, which includes a factor b 2 in the square of the wing area (the 
Sears 1 function also included is independent of b) , the mean-square 
total lift can be integrated in the form 




(b — ><») 


where Fj_, Fq 4 0. 

In summary, if the wing span is increased from zero while keeping 
all other conditions unchanged, at very small spans the mean-square lift 
first increases proportionally to the square of the wing area. For large 
values of the span, the low-frequency part of the lift spectrum approaches 
asymptotically a constant nonzero value, while the high-frequency part 
becomes asymptotically proportional to the span. Correspondingly, the 
asymptotic expression for the total mean-square lift has a constant term 
and a term proportional to the span. 


Discussion of Results 

The experimental results obtained in the present investigation are 
plotted and compared with the corresponding theoretical predictions in. 
figure 24. As pointed out previously, the turbulence scale A appearing 
in the derivation of the span correction factor g should be determined 
on the basis of the measured one -wave -number spectrum. Attempts were 
made to do so; however, since the quantitative values of the power spec- 
trum were doubtful and in any case extended over a limited frequency 
range, a reliable estimate of A could not be obtained in this manner. 
Dryden, Schubauer, Mock, and Skramstad (ref. 23) have investigated the 
turbulence scale downstream of coarse grids . If one neglects differences 
in the shape of grid members and possible effects of the background turbu- 
lence in the wind tunnel, the measurements of reference 23 would indicate 
for the present case a scale roughly 20 to 25 percent of the grid mesh ^ 
size, that is, A = 1.2 to 1.5 inches. This would mean 0.4 < A/b < 0-5, 
which is in good agreement with the experimental results at reduced fre- 
quencies above approximately 0.5; at lower frequencies the scatter pre- 
vents definite conclusions. 


Concurrently with the present investigation, Lamson (ref.. 13) has 
measured the aerodynamic transfer function of an essentially similar 
configuration. His results cover two turbulence scale-to-span ratios 
(10/7 and 5/7) and also demonstrate the effect of end plates shielding 
the measuring element from turbulence along the fixed span. The corre- 
sponding theoretical transfer functions, based on the present analysis, 
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are compared with Lamson's data in figure 28(a) (no end plates) and 28(h) 
(with end plates). In the former case, agreement in order of magnitude 
and dependence on turbulence scale-to-span ratio is obtained at k > 1. 

In case of end plates, the analysis could actually be expected to be more 
accurate in view of the strip-theory assumption; in fact, for k > 0.8, 
Lamson's measurements are now above the theoretical estimates. For lower 
values of the reduced frequency, the experimental results in both cases 
fall below the analytical curves and for k < 0.2 the dependence on tur- 
bulence scale-to-span ratio is actually reversed. 

Thus, at high-enough reduced frequencies, the present analysis pre- 
dicts an order of magnitude and dependence on turbulence scale-to-span 
ratio in essential agreement with existing experimental evidence . The 
quantitative difference between the measurements of Lamson and those of 
the present report may be traceable to factors such as the use of a wing 
of total aspect ratio 2.8 in an open jet (Lamson) and one of aspect 
ratio 5 in a closed wind tunnel (present investigation). Also, the 
method employed by Lamson for determination of turbulence scale is dif- 
ferent and actually makes use of measured characteristics of the power 
spectrum. 

It is believed that a more extensive, systematic series of measure- 
ments involving several wing configurations is still needed for complete 
and conclusive evaluation of the validity of the theoretical approach. 
Such measurements should be accompanied by a detailed investigation of 
the turbulence in order to determine its scale and to see how well the 
assumption of the approximate spectrum formula is justified. Analyti- 
cally, the weakest link appears to be the use of strip theory in the 
computation of the aerodynamic influence functions; theoretical research 
on this problem would be very worthwhile. 


CONCLUDING REMARKS 


Instrumentation for measurement and analysis of downwash and lift 
fluctuations was developed and found to operate satisfactorily. 

Experimental values of Sears' gust function could not be obtained 
because of excessive scatter of the measurements. The order of magnitude 
and trend of the measured downwash and lift, however, agreed with theo- 
retical predictions. Conclusive measurements would require improvement 
in environmental conditions, especially in the background turbulence 
level of the wind tunnel. 

Measurements were made of the transfer function between the power 
spectra of random downwash and the lift produced thereby on a simple 
rigid airfoil. An approximate theoretical result was also derived an 
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compared with the present measurements as well as with those of another 
investigation. Reasonable agreement was obtained at high reduced fre- 
quencies. Further experiments are recommended in order to evaluate com- 
pletely the validity of the theoretical approach. Improvements in the 
theoretical treatment would also be desirable , especially in regard to 
the aerodynamic influence functions . 


Massachusetts Institute of Technology, 
Cambridge, Mass., July 1, 1955* 
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APPENDIX A 

OPERATION OF x-WIRE HOT-WIRE PROBE FOR CROSS -VELOCITY 

MEASUREMENT 


The customary practical representation of the response character- 
istics of an oblique hot-wire is based on Kings 1 law combined with the 
assumption that the wire responds only to the velocity component normal 
to it. (See sketch and, for instance, ref. 25 *) 


v 

u U 


The resulting mathematical expression is 




l2 = 1 2 

x o 


2 a w' 

1 + V 


|U cos 


U r 


+ 1 


(Al) 


where I Q and U Q are calibration constants, I is the heating current, 
U is the velocity, \|r is the sweep angle, and a^' = (R w - the 

overheat factor or the relative resistance increase due to heating of the 
wire. In terms of the current I and cold resistance R e , the voltage 
output v of the wire can be expressed 


v 



(A 2 ) 


If the velocity is perturbed slightly by disturbance components u 
and v, the linearized perturbation in v will be 
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Av = 


ip 2 R e 

I 2 U 


1 + 


I 


U cos \|r 


1 + 


JJo 

U cos \|r 


i 1 

i_ 


I 2 


1 + 


l U cos 


-fu + -a—) 

>\ cot ty 


Cu + Dv 


(A3) 


provided that current I is kept constant, which is the case in the 
present hot-wire circuit. 

In an x -probe, another wire is placed at an angle of approximately 
-i|r with respect to the flow . If one now combines the outputs of the 
two wires as follows, applying a voltage stepdown device with factor r 
to one wire, the total response is 


AVi 2 = Av-j^ - Av 2 


= ^C}U + D^v^ - 
= - rC 2 )u + 


r ^C 2 u - D 2 vj 
+ rD 2 ^v 


(A4) 


If the wires and their operating conditions were exactly identical, one 
would have Cj_ = C 2 , and with no voltage stepdown (r = 1) the probe would 
respond to v-fluctuation only. 


Generally, it is not possible to produce two wires with exactly 
identical calibration constants. However, even in that case, the fac- 
tor C x - rC2 can be made equal to zero independently of the mean 

velocity U by suitable choice of the heating currents, the probe 
orientation with respect to the mean flow, and the stepdown ratio r. 
By examination of equations (A3) and (A4), one can write the following 
conditions : 


Heating-current ratio, 

ii = Mi 
^2 ( ^o) 2 


(A5) 


stepdown factor , 

( I °)l( R e)l 
( ■*■0) 2 ( R e) 2 


(a6) 


r 
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and probe orientation, 


cos _ (^°)l 
cos + 2 (U Q ) 2 


(AT) 


where ^ 2 . + ^2 i s cons ’* :an ’t* 

The calibration constants Iq and can be readily obtained in 

the form of equation (Al) from static calibrations of the individual 
wires, thus making possible the fulfillment of the first two conditions. 
The third one implies knowledge of the total angle between the two hot- 
wires. There is, however, no need to perform elaborate precision meas- 
urements and adjustments of the probe because the fulfillment of the 
three conditions (eqs. (A5) to (AT)) also make the total direct -current 
output 2 of ^e two-wire combination equal to zero, as can easily 
be seen in* the following equation: 


^ 1,2 = V 1 " rv 2 = 


MMi 


rlafe) 


1 - 


*T. , 

2(lo)i 2 ' 


U cos ^ 

( U o)l 


+ 1 


1 - 


2 fio) : 


U cos ^2 


( u °) 


+ 1 


(A8) 


Thus, one need only to adjust the heating currents and the stepdown 
voltage divider to satisfy equations (A5) &nd (A6); thereafter, the third 
condition (eq. (A7)) can be satisfied by turning the probe relative to the 
main velocity direction (i.e., by varying iq and \|r 2 ) such a posi- 

tion is found where the total direct -current output (eq. (A8)) vanishes. 

This procedure was successfully followed in practice. After adjust- 
ments of the probe were completed, its sensitivity to streamwise 
u-fluctuations was checked by varying the tunnel speed. Generally, the 
response to a u-disturbance was found to be at most a few percent of the 
response to a v-disturbance of the same magnitude. 
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APPENDIX B 


STATISTICAL RELIABILITY OF POWER-SPECTRUM MEASUREMENTS 


Since it is impossible to measure in practice the infinite time 
average implied by the use of generalized harmonic analysis, determina- 
tion of sufficient observation time to obtain values of power spectral 
density with a specified statistical reliability poses an important prob- 
lem. In general, the observation time has to be large compared with any 
characteristic times of the random process. However, since these quan- 
tities are themselves objects of measurement, being involved as param- 
eters in the power spectrum or correlation functions, one has to resort 
to assumptions regarding the probable shape of the power spectrum. 

Tukey (ref. 2 6) considered the accuracy of power spectra computed 
from a discrete set of sample points of the random function, including 
different spectrum shapes . The present analysis deals with the relia- 
bility of continuous mean-square measurements of the output of a narrow- 
band-pass filter subjected to the random input and is based on the line 
of thought presented by Rice (ref. 7). Davenport, Johnson, and Middleton 
(ref. 27) discuss the same problem employing assumptions somewhat dif- 
ferent from those of the present investigation. 

It is assumed that (l) the power spectral density to be measured is 
constant over the filter range; (2) the random process is Gaussian; 

(3) the band-pass filter is approximately a single tuned circuit with 
response 


( 1 +) -the squaring device is a heater-thermocouple unit with time lag; and 
(5) the averaging is done by direct integration of the squared signal. 

Let f(t) represent the filter output; then, the thermocouple out- 
put is 



t'-t 



,t 


e 
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where T c is the time constant of the thermocouple. The infinite time 
average of s is s = f 1 2 (ref. 7), while measurement over time T 
yields 


1 ■ * lo S(t)dt = * fo &{[1 


t'-t 


e Tc f 2 (t')dtjdt 

It is desired to determined the variance of s , that is, the quantity 




Now, 


T 


- s\2 _ 's‘ 2 - 

= P 


— r>T 


~ 2 = ±]j s(ti)s(t 2 )dt 1 dt 2 = J o (T - T)cp ss (T)dT 

where cp ss is the normalized correlation function of s(t). The prod- 
uct s 2 cp (t) can be derived from the filter response; by Fourier trans- 

form one finds that the normalized correlation function of the filter 
output is 


cpff(T) = e 

and according to Rice (ref. 7 ) 


2Q 


COS OCIqT 


s^P s (t) = s 


-2 


, 1 

t 

r°° " t c 

1 + — 


Tc, 

Jo 


cp ff 2 (t + t) + <p ff 2 (t - t)| dt 


The resulting mathematical expressions become rather lengthy; however, 
a convenient simplification can be made by assuming that Q 2 » 1 a^id 

(o^)T c ^ 2 » 1 which is true in many practical cases. By carrying out 
the detailed computations and successive substitutions one finally 
obtains (denoting o^T/Q = p and o^Tc/Q = 7) 


1 J 2y(l - 7 + 7 e~ p /^ - e~ p ) 

" P J 1 _ r 2 
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In the present investigation Q = 80, T c = 5 seconds, and 

3 < ^ < 150 cps; hence, approximately, 

2rt 


1 < 7 < 50 

Curves of cr(p) are plotted in figure 25 for 7=1, 10, and 50. It is 
seen that, for example, statistical accuracy of 3 percent root mean square 
implies the use of integration times of the order p = 50, that is, from 
5 seconds at 2it = 150 cps to over 3 minutes at <%/ 2jt = 3 cps . In 

practice, these times were exceeded by a substantial margin to make some 
allowance for deviations from the assumptions employed in the analysis. 

It is to be emphasized, however, that the 3 percent is the root-mean- 
square error of an infinite ensemble of measurements and, therefore, it 
is impossible to guarantee the accuracy of an individual measurement 
within specified limits. 

Some comment may be in order regarding the curves presented in fig- 
ure 25. At a low value of p the error decreases with increasing time 
lag 7, which is consistent with physical intuition about the effects 
of smoothing; however, at higher values of p the trend is reversed. 

The reason for this can be seen by considering the error equation, which 
may be written 


° 2 - £$)£ (t ■ T)<p ^ (T>aT - 1 

increasing time lag decreases s^/ s^ it also makes the fluctuations 
more correlated and thus increases 9 ss (t) in the range of integra- 

At long integration times, the latter effect becomes dominating. 

The values of a at p = 0 represent the standard deviation of 
an instantaneous observation and thus can be employed to estimate the 
frequency above which integration is no longer necessary . 
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Figure 1.- Side view of 5- by T.^-foot wind tunnel. 
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Figure 2.- Hot-wire probe. 
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Figure 3.- Hot-wire feeding and bridge circuits . Galvanometer selector, 
(l) - (l)' to (D - (D'; amplifier selector, (y - © to (8) - (g; . 


MCA TN 3878 




NACA TTJ 5878 


39 



Figure 4.- Typical single hot-wire calibration. 
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Figure 5 *- 


Typical cross -velocity calibration of double-wire probe. 
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Figure 7 .- Circuit diagram of strain-gage force pickups. 
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Figure 8.- Static calibration of strain-gage force pickup. Average sen- 
sitivity, 2.63 mv/v-kg; maximum deviation, 2.7 percent. 
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Figure 9.- Dynamic response of lift-pickup system to sinusoidal motion 
of support at fixed acceleration amplitude . 
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Figure 10.- Amplifier circuit (one channel). 
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Figure 11 .- Calibration of amplifier (push-pull, maximum gain). 


NACA TN 3878 


47 



3an±ndi/\iv indNi 


aarundwv indino 


Figure 12 . - Frequency response of amplifiers . 
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Figure 13.- Circuit for thermocouple driver. 
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Figure 14 . - Thermocouple calibration. 
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Figure 15. - Integrator calibration. Average slope, O.556 v/mv-min. 
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Figure l6.- Oscillator wing. 
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Figure 17.- Block diagram for measurements of sinusoidal downwash. 



(a) Complete signal. 



(b) Turbulence and amplifier noise. 

Figure 18.- Sample recordings of downwash signal. Frequency, 4.24 cps; 
recording sensitivity, 70 mm/deg flow deflection. 
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Figure 19.- Measurements of sinusoidal downwash. U = 50 mph- 
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Figure 20.- Block diagram for measurements of sinusoidal lift. 
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Figure 21.- Measurement of sinusoidal lift. 
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Figure 22.- Turbulence -generator grid. 
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Figure 23 .- System diagram for power-spectrum measurements . 
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Figure 2 b.- Comparison of experimental results and theoretical estimate 
of complete transfer function, h/c = 1/3. 
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Figure 25.- Statistical root-mean-square error 
ment with filter and time lag. 
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Figure 26 .- Coordinate systems and sketch of lift influence function, 
xy system stationary; |t] system follows wing moving with velocity 
in x-direction. 
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(a) With no end plates. 

Figure 28.- Comparison of results of reference 13 to theoretical estimate 

(b/c) = 0.28) . 
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